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We discuss the time dependence and future of the Cosmic Microwave Background (CMB) in 
the context of the standard cosmological model, in which we are now entering a state of endless 
accelerated expansion. The mean temperature will simply decrease until it reaches the effective 
temperature of the de Sitter vacuum, while the dipole will oscillate as the Sun orbits the Galaxy. 
Ifowever, the higher CMB multipoles have a richer phenomenology. The CMB anisotropy power 
spectrum will for the most part simply project to smaller scales, as the comoving distance to last 
scattering increases, and we derive a scaling relation that describes this behaviour. However, there 
will also be a dramatic increase in the integrated Sachs- Wolfe contribution at low multipoles. We 
also discuss the effects of tensor modes and optical depth due to Thomson scattering. We introduce 
a correlation function relating the sky maps at two times and the closely related power spectrum 
of the difference map. We compute the evolution both analytically and numerically, and present 
simulated future sky maps. 
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I. INTRODUCTION 

The Cosmic Microwave Background (CMB) radiation 
provides us with a vital link to the epoch before the 
formation of distinct structures, when fluctuations were 
still linear and carried in a very clean way information 
about their origin, presumably during a phase of infla- 
tion. The simple dynamics of the generation and prop- 
agation of CMB anisotropies (see e.g. Refs. [1, 2, 3] and 
references therein) depends on a handful of cosmological 
parameters, P^, such as the Hubble constant, the mat- 
ter density, and spatial curvature, in addition to the ini- 
tial conditions set through inflation. These dependencies 
have been thoroughly investigated over the past couple of 
decades and form the basis for estimating the parameters 
from the observed anisotropy spectrum of the CMB. 

However, there is one dimension in the parameter space 
of the CMB that has received little explicit attention. For 
flxed matter content and curvature of the Universe today, 
we still have the freedom to evolve the CMB anisotropies 
forwards or backwards in time. For the practical business 
of performing CMB parameter estimation, it is natural of 
course to suppress this freedom, since we are interested 
in predicting the anisotropies today. The constraint to 
"today" can be applied in at least two ways, which it 
is important to distinguish. From the set of parameters 
Pi we can calculate the proper-time age of the Universe, 
to- This quantity is only determined to an accuracy set 
by the parameters Pi (e.g. using Wilkinson Microwave 
Anisotropy Probe (WMAP) 3-year data, Spergel et al. 
[4] find that to = 13.73t|^;i? Gyr). However, the WMAP 
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results constrain the redshift of last scattering, Zioc, (de- 
fined as the centre of the recombination epoch) to much 
greater accuracy: z-^^c = 1088^2 [5]. This very small un- 
certainty is the result of our accurate determination of 
the mean temperature of the CMB, T = 2.725 ± 0.001 
K [6, 7]. Thus, even though to is only known to an ac- 
curacy comparable to the other parameters P^, implicit 
in analyses of the CMB is the very tight constraint on a 
different temporal coordinate, Zjoc or T. 

Essentially, the constraint on to arises from our deter- 
mination of the expansion rate today, together with in- 
formation on the content and geometry of the Universe, 
which affect its expansion history. The constraint on T 
is entirely independent of the content or geometry, hence 
its superior accuracy. Popular CMB anisotropy numeri- 
cal packages, such as CMBFAST or CAMB [8, 9] [38] impose 
the tight constraint arising from the mean temperature 
through an input parameter. This constraint on T is 
equivalent, via the Stefan-Boltzmann law, to a constraint 
on the energy density in CMB radiation, p-^. Therefore 
by adjusting the temperature input parameter, and other 
cosmological input parameters accordingly, it is possible 
to generate spectra with these packages that correspond 
to a given model evolved into the past or future. Alterna- 
tively it is possible to modify the code in these packages 
to directly integrate to future times without the need to 
modify the input parameters. Note that if we just vary 
the proper age to of a model by varying the expansion 
rate today, this necessarily changes the relative contribu- 
tions of matter and radiation in the past, which affects 
the physics at recombination and hence the shape of the 
CMB spectrum. 

Fortunately the required code modifications are rela- 
tively straightforward, and in addition it is possible to 
describe the temporal evolution of the CMB anisotropies 
analytically to very high precision. In this work we sys- 
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tcmatically describe this evolution both numerically and 
analytically, within the context of the standard ACDM 
(A Cold Dark Matter) model, in order to complete the 
standard results on the parametric dependence of the 
CMB. The verification of our numerical work with our 
analytical results, and conversely the characterization 
of our analytical approximations with the full numeri- 
cal calculations, will be crucial in this novel study. We 
will find that while the temporal behaviour of the CMB 
power spectrum is determined mainly by a simple ge- 
ometrical scaling relationship, less trivial physics arises 
when we consider the behaviour of correlations between 
anisotropics at different times. 

It can certainly be argued that the standard calcu- 
lations of the CMB anisotropy spectrum implicity de- 
scribe its time dependence in that the spectrum must be 
evolved from the time of recombination to the present. 
Nevertheless, there appear to have been very few explicit 
discussions of the time dependence, with the exceptions 
being primarily concerned with the distant future. Gott 
[10] points out that at extremely late times the typical 
wavelengths of CMB radiation will exceed the Hubble 
radius, and so the CMB radiation will be lost in a de 
Sitter background. Loeb [11] mentions that as the time 
of observation increases, the radius of the last scattering 
sphere also increases, and approaches a maximum in a 
ACDM model. This leads to the potential for reducing 
the cosmic variance limitation on the determination of 
the anisotropy spectrum. Krauss and Scherrer [12] point 
out that well before this final stage, the CMB will redshift 
below the plasma frequency of the interstellar medium 
and hence be screened from view inside galaxies. 

Importantly, when discussing the distant future evo- 
lution of the Universe it must be remembered that even 
the qualitative details can depend very sensitively on the 
model adopted. An extreme example is the potential de- 
struction of the Universe in finite proper time in a "big 
rip" [13], when the dark energy violates the weak energy 
condition, with equation of state w < —1. In the present 
work, for the sake of definiteness and simplicity, we con- 
servatively choose a spatially flat model in which the dark 
energy is a pure cosmological constant, with cosmological 
parameters consistent with the WMAP results [4] . How- 
ever, using the techniques we discuss it is straightforward 
to extend our results to other specific dark energy mod- 
els. 

An interesting question that naturally arises in the 
present context is: How long must we wait before we 
could observe a change in the CMB? The formalism that 
we develop here will be necessary to answer this question, 
and we will address this explicit issue in separate work 
[14]. We stress that the practicalities of an experimental 
search for the effects we describe is not our concern in 
the present work. In addition to setting the stage for the 
detectability question, we believe that a discussion of the 
evolution of the CMB sky is useful in its own right in 
elucidating the physics of the anisotropics in a vacuum- 
dominated cosmology. 



We begin in Section II with a description of the time 
dependence of the "bulk" properties of the CMB, namely 
the mean temperature and the dipole. After a brief re- 
view of the formalism used to describe the anisotropy 
power spectrum, its evolution is described analytically in 
Section HI, including the effects of the integrated Sachs- 
Wolfe effect, tensor modes, and reionization, and numer- 
ical calculations arc presented using our modified version 
of the line-of-sight Boltzmann code CAMB. In Section IV 
we introduce the difference map power spectrum and as- 
sociated correlation function, and present analytical and 
numerical calculations. Section V presents our conclu- 
sions, and in the Appendix a description of an impor- 
tant approximation method is presented. We set c = 1 
throughout. 

II. TIME EVOLUTION OF THE BULK CMB 

The temperature fluctuations on the CMB sky can be 
decomposed into a set of amplitudes of spherical harmon- 
ics (see Section HI A). The angular mean temperature 
(or "monopole") and dipole have a special status. The 
mean temperature is just a measure of the local radia- 
tion energy density, while the value of the dipole depends 
on the observer's reference frame at linear order (higher 
multipolcs are independent of frame at this order). 

A. The mean temperature 

As time passes, the change in the CMB that is simplest 
to quantify is the cooling of its mean temperature T due 
to the Universe's expansion. The CMB radiation was 
released from the matter at the time of last scattering, 
when T ~ 3000 K. It later reached a comfortable 300 K 
at an age of about < ~ 15 Myr, and is now only a frigid 
few Kelvin. Indeed, the monotonicity of the function 
T{t) means that T itself can be used as a good time 
variable. Thus wc can consider measurements of T as 
direct readings of a sort of "cosmic clock" . 

Today the CMB radiation is essentially free stream- 
ing, i.e. non-interacting with the other components of 
the Universe. Therefore the energy density in the CMB 
evolves according to the energy conservation equation 
p-y = —AHpy, where H is the Hubble parameter and 
the overdot represents the proper time derivative. Since 

(X T'^, we have T = ~HT. Evaluating this expression 
today, using Tq = 2.725 K and Hq = 73 kms-^Mpc"! 
(subscript o indicates values today), we find 

to = -0.20 pK kyr^\ (1) 

Thus in 5000 yr the mean temperature will drop by 1 pK. 

The CMB radiation continues to redshift indefinitely 
as the Universe expands in the late A-dominated de Sit- 
ter phase. However, this does not mean that as the CMB 
becomes increasingly difficult to measure, clever experi- 
mentalists need only to ever refine their instruments in 



3 



order to keep up. Instead, a fundamental limit exists 
below which a CMB temperature cannot be sensibly de- 
fined. An object in an otherwise empty de Sitter phase 
will see a thermal field with temperature [15] 



H 



1 



(2) 



Therefore, after the CMB temperature redshifts to below 
Tds, the CMB becomes lost in the thermal noise of the 
de Sitter background, as pointed out in [10] (see also 
[16]). 

To see explicitly the difficulty with measuring the CMB 
at such extremely late times, consider the typical wave- 
lengths of radiation in the CMB. A thermal spectrum at 
temperature T consists of wavelengths A on the char- 
acteristic scale (the precise peak position of the 
Planck spectrum depends on the measure used for the 
thermal distribution). Therefore when T = T^s wc have 
A ~ H~^, i.e. the typical CMB wavelengths become of 
order the Hubble length. Alternatively, at late times in 
the de Sitter phase the frequency of a mode of CMB ra- 
diation of fixed comoving wavelength redshifts according 
to 



w(0 =c^(ti)e-^^^(*-*i\ 



where 



^^dS = V ^aHq 



(3) 



(4) 



is the asymptotic value of the Hubble parameter and ti is 
some late proper time. Therefore the accumulated phase 
shift between time ti and the infinite future is 



uj(t)dt 



(5) 



This expression tells us that when the frequency becomes 
less than the Hubble parameter (i.e. the wavelength be- 
comes larger than the Hubble length), a full temporal 
oscillation cannot be observed, even if we observe into 
the infinite future. In terms of conformal time, the oscil- 
lation rate remains constant in the de Sitter phase, but 
there is only a finite amount of conformal time available 
in the future. Indeed, considering the quantum nature 
of such a mode, this calculation provides insight into the 
necessity of a residual de Sitter thermal spectrum at this 
scale. 

The energy density in the CMB at arbitrary scale fac- 
tor a is given by 



where mp is the Planck mass and fl^ ~ 5 x 10" 
fraction of the total density in the CMB today. Using 
this expression we can show that we must wait until 
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a/ao ~ 10^° before p-y 
comes lost in the de Sitter background. This corresponds 



PdS 



T^g and the CMB be- 



to an age oi t = 1 Tyr. If we ask instead at what scale 
factor would the radiation density be equal to the Planck 
density, pp = mp, we find a/ao ^ 10^'^^. It might ap- 
pear, therefore, that we exist at a special time, in that the 
radiation density today is roughly 120 decades removed 
from both the Planck era and the final era when T ~ T^s- 
To understand the origin of this coincidence, note that by 
virtue of the above expressions and the energy constraint 
(or Friedmann) equation, the three densities pp, pA, and 
Pds are in the geometrical ratio nip : Arrip : A^ , up to nu- 
merical factors. Therefore, the apparent coincidence just 
described is actually equivalent to the standard coinci- 
dence problem, namely that pA — Ptot today, given that 
Pj differs from ptot today by "only" a few decades. Any 
density that today even crudely approximates the dark 
energy density will necessarily be separated by roughly 
120 decades from both pp and pds- 



B. The dipole 

The observed dipole anisotropy in the CMB can be 
attributed to the Doppler effect arising from our peculiar 
velocity, v, with respect to the frame in which the CMB 
dipole vanishes. That peculiar velocity, and hence the 
dipole, is expected to evolve with time. The magnitude 
of the dipole can be specified by the maximum CMB 
temperature shift over the sky, ST^, due to the velocity 
V. This is given by the lowest order Doppler expression. 



STd 
T 



V. 



(7) 



(In terms of the spherical harmonic expansion to be intro- 
duced in Eq. (15), we have 6Td/T = ^3/(47r)aio, when 
the polar axis is aligned with v.) 

The current best estimate of the magnitude of the 
dipole comes from observations of the WMAP satellite — 
indeed, the annual modulation by the Earth's motion 
around the Sun is actually used to calibrate satellite ex- 
periments, so this aspect of the time- varying dipole is al- 
ready well determined. The measured value of the dipole, 
in Galactic polar coordinates, is {STd, I, b) = (3.358 ± 
0.0017 mK, 263.86 ± 0.04°, 48.24 ± 0.10°) [17]. Equiv- 
alently, the Cartesian velocity vector is vq ~ (—26.3, 
—244.6, 275.6) kms^^, where the first component is to- 
wards the Galactic centre, and the third component is 
normal to the Galactic plane. Therefore, in natural units 
we have vq — 1.2 x 10"'^, so the lowest order approxima- 
tion, Eq. (7), is valid. 

In order to determine the evolution of the dipole, we 
could straightforwardly calculate it at linear order. How- 
ever, linear theory is certainly not a good approximation 
on sufficiently small scales today. To fully describe the 
evolution of the velocity v we must take into account the 
presence of the nonlinear structures we observe on small 
scales today. 

This velocity vector can be considered as a sum of in- 
dividual vectors contributing to the overall motion of the 
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Sun with respect to the CMB. In the local neighbour- 
hood, the Sun moves with respect to the "local standard 
of rest" , which in turn moves with respect to the Galac- 
tic centre. However, the peculiar motions in the Solar 
neighbourhood are at the 10% level compared with the 
motion of the Sun around the Milky Way [18], so for 
the purposes of the simple calculation which follows we 
ignore these contributions. Also, we will consider here 
time scales short enough that the motion of the Milky 
Way within the Local Group, and the Local Group rel- 
ative to Virgo, the Great Attractor, and other distant 
cosmic structures is approximately constant (see, e.g., 
[19] for a description of these motions). 

Just as today we can detect the modulation of the 
Earth's motion around the Sun, in the future, with in- 
creasing satellite sensitivity, we may be able to observe 
the Sun's motion around the Galaxy. For the motion of 
the Sun around the Milky Way, we assume that this is 
simply a tangential speed of 220 kms~^ at a distance of 
8.5 kpc. Using the current observed value of v to infer 
the velocity of the Galactic centre with respect to the 
CMB rest frame, the time dependent Sun-CMB velocity 
vector is then 



vit) 



/2nt\ 
222 sin — - 26.3, 



222 cos 



466.6, 275.6 



kms^S (8) 



where the Galactic orbital period is T = 2.35 x 10^ yr. 

In order to ascertain when a change in the dipole is 
detectable, one could compute a sky map of the dipole at 
two times. If the temperature variance of the difference 
map is greater than the experimental noise variance, then 
a detection is probable. In this case, the variance of the 
difference map, which we denote by Cs, is Cs = [{5vxY + 
[SvyY + ('5'^z)^]/(47r), where 5v is the difference of the 
Sun-CMB dipole vector between the two observations. 
Using Eq. (8), and converting to fractional temperature 
variations, the signal variance of the changing dipole is 
then 



Cs 



8.7 X 10"^ 



2lTt 



(9) 



Later in this paper we compute signal variances involving 
higher order CMB multipoles. These variances are of 
course much smaller than that of the dipole. In a follow 
up paper [14] we will discuss in detail the prospects for 
detecting a change in the CMB with future experiments. 

Finally, we note that in this simple calculation it is 
assumed that we have a frame of reference external to the 
Milky Way in order to construct our coordinate system. 
This could be provided, for example, by the International 
Celestial Reference Frame, based on the positions of 212 
extra-galactic sources [20]. 



III. THE ANISOTROPY POWER SPECTRUM 
A. Review of the basic formalism 

There is much more information encoded in the 
anisotropies of the CMB than in the mean temperature, 
since the anisotropies are determined by the details of the 
matter and metric fluctuations near the last-scattering 
surface (LSS) and all along our past light cone to to- 
day. Therefore it is much less trivial to determine the 
time evolution of the anisotropies than the mean tem- 
perature (or dipole). However, in the approximation that 
all of the CMB radiation was emitted from the LSS at 
some instant ^ls when electrons and photons decoupled, 
and then propagated freely, the evolution of the primary 
power spectrum of the CMB is determined by a simple 
geometrical scaling relation which is closely related to the 
main geometrical parameter degeneracy in CMB spectra. 
In order to derive this relation, and to describe the be- 
haviour of the correlation functions introduced in later 
sections, it will be helpful to first summarize the stan- 
dard description of CMB anisotropies in a form that will 
be easy to generalize. This subsection may be skipped 
by readers familiar with the material. For detailed treat- 
ments of the generation of anisotropies see e.g. [21, 22]. 

At very early times, when each perturbation mode, la- 
belled by comoving wavevector fc, is outside of the Hub- 
ble radius, the fluctuations can be described by a single 
perturbation function, for the case of adiabatic pertur- 
bations. It is very convenient to take this function to be 
the curvature perturbation on comoving hypersurfaces, 
72,, since this quantity is conserved on large scales in this 
case, and hence can be readily tied to the predictions of 
a specific inflationary model. In the simplest models of 
inflation, TZ is predicted to be a Gaussian random field to 
very good approximation, fully described by the relation 



{n*{k)nik')) = 2TT^s^{k-k') 



(10) 



with primordial power spectrum V-jzik) and k = \k\. For 
a scale-invariant spectrum we have V-ji{k) = constant. 

The fluctuations at last scattering can be described by 
a set of matter and metric perturbations, </>,;(a;, 77), where 
for future convenience we have used comoving coordinate 
X and conformal time r]. Since linear perturbation theory 
is a very good approximation at the scales sampled by the 
CMB, this set of perturbations is determined from the 
primordial comoving curvature perturbation by transfer 
fimctions Ai{k,'q) via 



0,(fc,77) A,(fc,77)72(fe). 



(11) 



In the approximation of abrupt recombination, so that 
the LSS has zero thickness, followed by free streaming of 
radiation, the observed primary temperature anisotropy 
5T{e)/T in direction e is determined by the fluctuations 
at the corresponding point on the LSS, i.e. 

5T{e,i^) 



^'(0j(?'LS,e,?/Ls)), 



(12) 
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for some linear function F. Here tls = V ~ VhS is the 
comoving radial coordinate to the LSS from the point 
of observation, taken to be the origin of spherical coor- 
dinates. Eq. (12) ignores both the effect of gravitational 
lensing by foreground structure and the effect of reioniza- 
tion at late times (we examine the effects of reionization 
in Section III E) . In the approximation that photons are 
tightly coupled to baryons before ryLs i the function F can 
be written in terms of two perturbation functions as 



^'(0i(''LS,e,77Ls)) = (/>i(?'LS,e,77Ls)- 



d 



drhs 



02(?'LS, e,77Ls)- 

(13) 

Eqs. (12) and (13) describe the generation of CMB 
anisotropics through the Sachs- Wolfe effect [23], with 
the first term on the right-hand side of (13) the so- 
called "monopolc" contribution, and the second term the 
"dipole" or "Doppler" contribution. 

The preceeding equations have the very simple inter- 
pretation that when we measure the CMB anisotropics at 
some time 77 we are "seeing" the primordial fluctuations 
TZ on the comoving spherical shell r = tls = V ^ '/ls 1 as 
processed by the linear transfer functions Ai and A2 to 
the time r/Ls- If wc observe at a later time 77', we see the 
fluctuations on a larger shell of radius r — r^g = rj' — ?7Lg, 
as illustrated in Fig. 1. The fluctuations at the LSS con- 
tain structure at various scales, encoded in the trans- 
fer functions, due to acoustic oscillations within the pre- 
recombination plasma. Assuming the statistical homo- 
geneity of space, that structure will occur at the same 
physical scales on the shells r = tls and r = r^g . There- 
fore we expect that structure visible at time 77 on angular 
scale 9 will also be visible at 77', but at the smaller angular 
scale 



rhs 



'LS 



(14) 



at least for small scale structure, 9-^1. To make this 
rigorous, and to derive in addition the scaling law for the 
amplitude of angular structure, wc need to next introduce 
a spherical expansion of the CMB anisotropy. 

We expand as usual the temperature fluctuation ob- 
served in some direction e in terms of spherical harmonics 

Yirn as 



Tiv) 



(15) 



im 



The expansion coefficients aim determine all the details 
of the particular sky map of the CMB observed at time 
77. However, the statistical properties of the aimS are 
determined through Eqs. (11) to (13) by the statistics of 
TZ encoded in Eq. (10). To make this explicit, we need the 
spherical expansion of the perturbations (/'i(7'LS, e, 7/ls)j 
namely 

0i(?'LS,e) = J kdk 4'iem{k)je{krLs)Yem{e). 

(16) 




FIG. 1: Our spheres of last scattering at time 77 (inner) and 
77' (outer). The group of vertical lines indicates the crests of 
a mode fe, and hot spots will be observed at the intersections 
of those crests with the spheres. One wavelength of the mode 
will span angle at 77 but a smaller angle d' at 77'. 



Here i = 1 or 2, is the spherical Bessel function of the 
first kind, and we have dropped the understood argument 
77LS- Next the identity 



f,m{k)=ki' / dnkf{k)Y,m{k) 



(17) 

(see, e.g., Ref. [24]) combined with Eq. (11) allows us to 
write 

(j)Um{k)^ A,{k)nira{k). (18) 

Now, combining Eqs. (13), (16), and (18), we have 

^(^LS,e) = / fcdfc^7^,™(fc)T(fc,^,rLs)i£™(e), 

(19) 

where 

r(fc,£,rLs) = Ai(fc)j,(fc7^s) +^2(fc)jKfcrLs) (20) 

and the prime denotes differentiation with resepect to 
tls- Finally, equating coefficients between Eq. (19) [with 
Eq. (12)] and Eq. (15), we obtain 



aha{.ri) = \ - kdknim{k)T{k,i,n^s), (21) 



where we have restored the argument 77 = t^ls + ??ls- 
This expression gives the CMB anisotropy in terms of 
the primordial perturbations TZ and a new linear transfer 
function T{k,£, tls)- 
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In order to determine the statistical properties of the 
aimS, we need the expression 

(7^,™(fc)7^,^„,(fc')> =27r25(fc-fc')^|P<5«'<5w, (22) 

which can be derived from Eq. (10). Using this expression 
and Eq. (21) we find 

{atmiv)ai'm'{v)) ^ Ci{ini)5u,5mm', (23) 

where 

^VTzik)T'(kJ,r^s)- (24) 

That is, each coefficient aim has variance Ci{rf) (which 
is called the anisotropy power spectrum) and coefficients 
for different spherical modes are uncorrelated. 

Note that Eq. (21), and hence Eqs. (23) and (24), hold 
even when we relax the tight coupling and free streaming 
approximations, with some transfer function T(k, tls). 
However, in the general case the transfer function must 
be calculated numerically. 

B. Analytical time evolution for primary 
anisotropies 

The formalism developed in the preceeding subsection 
can now be applied to describe the time evolution of the 
primary CMB anisotropy spectrum, under the abrupt re- 
combination and free streaming approximations. To de- 
termine the time evolution of C|(r?), Eq. (24) tells us that 
we only need to consider the behaviour of T^(fc, £, tls) as 
tls increases (note that we will often adopt the coordi- 
nate tls as an effective time coordinate). To do this, 
Eq. (20) tells us that we only require the behaviour of 

the products j|(fcrLs), (fcj'Ls), and (fcrLs)j£(fcrLs) as 
functions of tls- This can be done in the limit ^ ^ 1 us- 
ing asymptotic forms for the Bessel functions. For large 
I we can write [see Ref. [25], Eq. (9.3.3)] 

ji{x) = (x^ - x^f)-^/^[cos{e) + ©(l/f)], for x>e, 

(25) 

where = 6{x) is a real phase. For x < I, je{x) decays 
rapidly. This allows us to write a scaling relation for the 
envelope of the Bessel oscillations, namely 

ji{x) ^ ajai{ax), (26) 

for positive a such that 3> 1 also applies. This expres- 
sion will allow us to obtain the time dependence of the 
"monopole" contribution to T^(fc,£, tls), which is pro- 
portional to j|(fcrLs)- (The Bessel oscillations arc rapid 
relative to the range of scales that contribute to the inte- 
gral Eq. (24) for i :s> 1, and hence can be ignored.) We 
can write the "dipole" part j^^(fcr-Ls) in terms of spheri- 
cal Bessel functions using recurrence relations and again 
apply Eq. (26) to obtain the time dependence. The cross 



term proportional to je{kri^s)j'^{kri^s) can be shown to 
be negligible, i.e. the monopole and dipole contributions 
add incoherently. Applying Eq. (26), then, we find that 
for large £, 

r2(fc,£',rLs)^ %T2(fc,£,rLs), (27) 
where we have defined 

r' 

£' = i-^. (28) 

Applying Eq. (24) we finally obtain the scaling relation 
for the power spectrum, 

£'^a'{'n')^fCi{f]). (29) 

Importantly, Eq. (29) holds independently of the form 
of the functions Ai{k) and Vnik), so the result applies 
to the acoustic peak structure as well as to non-scale- 
invariant primordial spectra. 

This result confirms our previous prediction, Eq. (14), 
for the dependence of angular scales on observation time. 
But the dependence of the amplitude of the spectrum en- 
coded in Eq. (29) is also not surprizing, since the quantity 
£{£ + l)Ce is independent of £ in the Sachs- Wolfe plateau 
for a scale invariant spectrum, as is well known. But the 
height of that plateau, calculated using Ai — const and 
A2 ~ above, is independent of the observation time. 
(Indeed that height is, up to numerical factors, simply 
Vti- Recall that Ci is determined by the ratio ST/T. 
The absolute anisotropies ST exhibit the same expan- 
sion redshift as does the mean temperature T.) Hence as 
77 increases, the quantity £'^Ci(t]) must remain constant 
(up to corrections of order 1/^), which is precisely what 
Eq. (29) says. Of course, the result (29) is valid for the 
entire acoustic peak structure, not just the Sachs- Wolfe 
plateau. 

The result (29) is derived in the Appendix much more 
directly, without resorting to properties of Bessel func- 
tions, using the flat sky approximation. In that approach 
we consider anisotropies in a patch of sky small enough 
that it can be approximated as flat, and errors are again 
of order l/£. 

In addition to the main temperature anisotropies we 
have been considering here, there are also polarization 
spectra present in the CMB radiation. The polarization 
is sourced primarily near last scattering, so its spectra 
will also scale according to Eq. (29). A small part of 
the largest-scale polarization is sourced near rcionization, 
so we expect that that contribution will scale with the 
comoving radius to the rcionization redshift, rather than 
to the last scattering surface. 

Having found the scaling relation (29), we can next 
derive some simple consequences from it. First, we can 
write the total power in the anisotropy spectrum as 

Y,Ci{i^)^Y.^2£+l)Ci{f^)^2 f£Ce{v)d£ (30) 
em i '' 
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in the large f, approximation. Then, using Eq. (29), we 
have 

^(2£+l)Q(r;)c.^(2£+l)CK^'), (31) 



where the approximation comes from ignoring terms of 
order What this expression says is that the total 

power is constant in time, for the free streaming of pri- 
mary anisotropics. This result is equivalent to the "con- 
servation condition" stated in [26]. Implicit in this re- 
sult is the assumption of statistical homogeneity, so that 
no new anomalous power will be revealed at the largest 
scales as tls increases. As we will see in Section IIIC 
below, secondary anisotropies, in particular those gen- 
erated through the "integrated Sachs- Wolfe effect", are 
expected to grow dramatically at late times and hence 
the total power will not in fact be conserved. 

Another consequence of Eq. (29) follows from the na- 
ture of the asymptotic future in our ACDM model. Ob- 
servers in a universe with positive A have a future event 
horizon, i.e. the conformal time converges to a finite con- 
stant 77/ as proper time i — > 00. Therefore the angular 
scaling relation (28) tells us that as proper time (or scale 
factor) approaches infinity, the £ value for any particular 
feature in the spectrum, such as a peak position, will 
approach a finite maximum, i.e. features will approach 
a non-zero minimum angular size. (Geometrically, the 
LSS sphere approaches a maximum comoving radius, so 
features on it must approach a minimum size.) For our 
fiducial model we chose = 0.77, and so Eq. (28) gives 
for the limiting scaling relation 



(da) ^da 
(da)^irfa 



(32) 



For example, the first acoustic peak, which we observe 
to be at the position £0 = 221 today, will asymptote to 
If = 290 in the late dc Sitter phase. This asymptotic be- 
haviour is in marked contrast to that of a purely matter- 
dominated Einstein-de Sitter model. In the vanishing A 
case, the numerator in Eq. (32) diverges (no event hori- 
zon exists) and the structure in the Cg spectrum shifts 
to ever smaller scales. 

The geometrical scaling relation (29) is very closely 
related to the well-known geometrical parameter degen- 
eracy in the CMB anisotropy spectrum between spatial 
curvature and A [27, 28, 29]: If two cosmological models 
share the same primordial power spectrum Viz{k), the 
same physical baryon and CDM densities today, pb and 
Pc and finally the same angular diameter distance (Ia 
to the LSS, then they will exhibit essentially identical 
primary Ce spectra. The degeneracy can only be broken 
by secondary sources of anisotropy, such as the integrated 
Sachs- Wolfe effect, or by other cosmological observations. 

To understand the origin of this degeneracy and its re- 
lation to the preceeding discussion, recall that the energy 
density in the CMB today, p^, is fixed to very high ac- 
curacy by the measurement of the mean temperature, as 



we mentioned in the Introduction. Therefore if we con- 
sider models with identical values of the densities p^ and 
Pc today, then the densities of baryons, CDM, and pho- 
tons at last scattering are the same for all such models, 
since the densities scale in a well-defined manner (for ex- 
ample, P'j/pc oc ao/a). Therefore, given the same initial 
conditions in the form of Vnik), models that have com- 
mon values of Pb and pc today will have identical local 
physics at least to the time of recombination, when any 
spatial curvature or A will have negligible effect. Thus 
these models will produce identical primary anisotropies. 

If the models have different values of ^lK (spatial cur- 
vature) and A, then the dynamics, including the propaga- 
tion of CMB anisotropies, will differ significantly at late 
times as those components come to dominate. However, 
if the models share the same angular diameter distance, 
then their d spectra, which should be calculated using 
Eq. (24) with tls replaced by dA (at least for small scales 
where the effects of spatial curvature on the primordial 
spectra can be ignored), will be identical. Geometrically, 
models with identical Vn{k), pb, and pc share the same 
local physics to recombination, and hence the same phys- 
ical scales for acoustic wave structures (in particular the 
same sound horizon) . For models which additionally have 
identical d^, observers see the anisotropics generated on 
a spherical shell at the time of last scattering of identi- 
cal physical surface area (given by 47rc?^). Hence those 
physical acoustic scales are mapped to identical angu- 
lar scales in the sky for the different models. In short, 
the observed primary anisotropies in models with identi- 
cal Vizik), Pb, Pc and dA are produced under the same 
local physical conditions on a sphere of identical phys- 
ical size, and hence appear identical. The scaling rela- 
tion (29) describes how the observed anisotropies change 
if we hold the local physics at recombination (together 
with A and ^k) constant, but allow the time of observa- 
tion to vary, which amounts to simply varying the size of 
the sphere at last scattering that generates the observed 
anisotropies. The parameter degeneracy states that the 
same anisotropy spectrum can be produced even if A and 
VLk vary, as long as the size of the last scattering sphere 
is held constant. 

To close this discussion of the primary anisotropies, 
we introduce the power spectrum difference dCe{ri) = 
Ci{rj') — Ce{ri) between the spectra observed at two dif- 
ferent times. This is a measurable quantity which we 
might consider a candidate for detecting the evolution of 
the CMB. Given some spectrum C£{r]) at a single time r] 
we can readily calculate the difference SCe{r]) using the 
scaling relation (29). For small Srj = rj' — rj, we have 



(33) 



so that the change in the CMB power spectrum at fixed £ 
is proportional to drj. As we will see in the next Section, 
this behaviour differs from that of the power spectrum of 
the difference aemiv') ~ o-imiv)- Using Eq. (29) for the 
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time dependence, we can write 



' de 



(34) 



at first order in Srj/rji^s. Note that dCi can have either 
sign, and will equal zero whenever d{l'^Ci)/d(. = 0, as for 
example on a scale-invariant Sachs- Wolfe plateau or at 
an acoustic peak. 

Recall that the quantity Ci{ri) describes the relative 
anisotropics 5T /T , and hence is insensitive to the bulk 
expansion redshift. If we wish to consider instead the 
evolution of the absolute temperature anisotropics 
the relevant quantity to calculate is 



-2- 



5i] 



(35) 



iaH)- 

at lowest order in 5ri/{aH)^^, where we have used the 
relation T = —HT. Since in standard ACDM models 
77LS is a few times the comoving Hubble radius {aH)~^, 
the first term on the left-hand side of Eq. (35) is of the 
same order as the second term. That is, the expansion 
cooling effect is on the same order as the geometrical 
scaling effect, and so it will be important to distinguish 
the two processes. 



C. Integrated Sachs- Wolfe effect 

The simple scaling relation derived in the previous sub- 
section determines the time evolution of the power spec- 
trum of anisotropics produced near the LSS. However, 
in the standard ACDM model, significant anisotropics 
arc also produced at late times as a result of the chang- 
ing equation of state as the Universe becomes cosmolog- 
ical constant dominated. This process is known as the 
(late) integrated Sachs- Wolfe (ISW) effect [23, 30]. Since 
these anisotropies are produced relatively locally, their 
time dependence must be explicitly calculated. Note 
that anisotropies are also generated by the early ISW ef- 
fect, during the time that radiation still significantly con- 
tributes to the dynamics. However, those anisotropies are 
produced relatively close to the LSS, adding coherently 
to the primary Sachs- Wolfe contribution, and hence scale 
as do the primary anisotropies, according to Eq. (29). 

The contribution of the late ISW effect can be de- 
scribed by adding to the transfer function, Eq. (20), a 
term Tisw{k,i,ri) which is an integral over the fine of 
sight to the LSS. For the case of interest, for which the 
anisotropic stress is negligible, we have [22] 

Tisw(fc,^,?7)=2A3(fc) r dr^'g'{v')n[k{r^-r,')]. (36) 

J f)LS 

Here g'{ri) = dg/drj, with g(r]) being the growth func- 
tion which describes the temporal evolution of the zero- 
shear or longitudinal gauge curvature perturbation, ip 
(also called the "Newtonian potential"), via 



The function A3{k) is defined such that g{r]) — > 1 at early 
times in matter domination. Then, a gauge transforma- 
tion between the comoving, TZ, and zero-shear, rp, curva- 
ture perturbations during the matter dominated period 
gives A3{k) = —3/5. 

To evaluate the ISW contribution, we need to first de- 
termine the evolution of the curvature perturbation -0. 
To do this, we only need to solve the space-space, or dy- 
namical, linearized Einstein equation. For the case where 
pressure and anisotropic stress perturbations can be ig- 
nored, which holds in a universe containing only dust and 
A, this equation becomes 



4) + AHil' + (3i?2 + 2H)'iJj = 



(38) 



(see, e.g., Ref. [31]). There are no spatial gradients in 
this equation, which confirms that the growth function is 
independent of fc. It is straightforward to verify that the 
growing mode solution to this equation is 



H 



(39) 



Next, using the relation a^H — const, which holds ex- 
actly for a universe consisting of dust and cosmolog- 
ical constant, employing Eq. (37), and matching the 
growing mode solution to the initial condition ipo{k) = 
— (3/5)7?.(fc), we obtain 



5 noH 



di]' 



(40) 



where fio is the density in matter today relative to the 
critical density, and the scale factor is normalized to unity 
today. 

Now that the evolution of the growth function has been 
determined, we can evaluate the ISW contribution to the 
power spectrum. It can be shown that the cross term 
between the Sachs- Wolfe and (late) ISW terms is negli- 
gible, so the two add incoherently in Cf - For the ISW 
part we have 

/dk 
j^Vniml^ikJ,r^) (41) 

''\iig'\v')iv-v')- (42) 



25 P 



^ik,r,)=g{r^)A3{k)n{k). 



(37) 



To obtain this result wc have assumed a scale invariant 
primordial spectrum, V-jiik) = V-ji, and we have used the 
relation [32] 

(43) 

which holds for large £. Unfortunately it is at small £ 
that the ISW effect is greatest, so this approximation, 
which appears to be the best we can do analytically, is 
not terribly accurate at very late times, when, as we shall 
see, the ISW contribution becomes very large. Neverthe- 
less, Eq. (42) will give a reasonable estimate of the ISW 
contribution to a change SCi over short time intervals. 
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Given a primordial amplitude, Pk, and a matter den- 
sity parameter today, Hq, Eqs. (40) and (42) allow us to 
calculate the ISW contribution to the power spectrum 
at any time. In addition, taking the time derivative of 
Eq. (42), we find for the rate of change of the ISW con- 
tribution 



drj ^ 



iv) 



72fVn 
25 P 



di^'g'^il'). (44) 



Therefore, combining this expression with Eq. (34) for 
the change in the primary Sachs- Wolfe power spectrum 
over short time intervals (5r/, we find for the total contri- 
bution 



St] 
■rjLS 

72 tt^Vtz _ 
25 £3 



de 



-r]LS 



2Ceir]) 



dv'g'Hv') 



(45) 



D. Gravitational waves 

Inflationary models generically predict a spectrum of 
primordial gravitational waves, although the relative 
contribution of these tensor modes to the total CMB 
anisotropy ranges from substantial to very small, depend- 
ing on the model. The anisotropics arise through the 
tensor analogue of the scalar ISW line of sight integral, 
Eq. (36). In place of the time derivative of the scalar cur- 
vature perturbation ip, the tensor contribution involves 
the rate of change of the transverse and traceless part 
of the spatial metric perturbation, hij. The evolution 
of this part of the metric perturbation is given by the 
dynamical Einstein equation, which becomes (see, e.g., 
[31]) 



(46) 



in the absence of tensor anisotropic stress, which is a valid 
approximation in the matter- or A-dominated regimes. 

The dynamics of hij as dictated by this equation de- 
pends on the mode wavelength relative to the Hub- 
ble scale. For k/{aH) ^ 1, the tensor mode is over- 
damped, and the (growing) mode decays very slowly. For 
k/{aH) ^ 1, the mode undergoes underdamped oscilla- 
tions, decaying like hij oc e*'"'''/a in the adiabatic regime. 
Therefore, for a particular tensor mode fc, the rate of 
change h[j is peaked near the time that the mode crosses 
the Hubble radius, k/{aH) ^ 1, and is small at early 
and late times. This means that the contributions to the 
CMB anisotropics from a particular scale arise primarily 
when that scale enters the Hubble radius. In particular, 
scales that enter significantly before last scattering will 
have decayed before they could source the CMB. This 
imposes a small scale cut-off on the tensor anisotropy 
spectrum, with negligible power for 



where 1 /cils is the comoving Hubble radius at last scat- 
tering. Since ols is fixed for a particular model, the ^ 
cut-off scales with time of observation according to 



(48) 



which is the same scaling as for primary features in the 
scalar CMB spectrum, Eq. (28). 

For £ <^ £c, detailed calculations for a matter domi- 
nated universe [33, 34] show that the tensor anisotropy 
spectrum is nearly flat, mimicing the Sachs-Wolfe 
plateau. In this case the tensor spectrum does not evolve 
apart from the scaling of the cut-off, Eq. (48). How- 
ever, as described above, the largest angular scales will be 
sourced at the latest times, and so for a universe with cos- 
mological constant this will lead to some dependence on 
observation time for those scales as the equation of state 
changes at late times. For example, for time of observa- 
tion 7] the tensor quadrupole is sourced near very roughly 
the conformal time 77/2 [33]. In the matter-dominated 
era, the comoving Hubble radius \/{aH) = \/a increases 
with time, but as the Universe enters the A-dominated 
era, l/{aH) starts to decrease (indeed this deflnes accel- 
eration). Therefore, the largest scale modes that con- 
tribute to the tensor anisotropics will enter the Hubble 
radius somewhat later in the presence of a cosmologi- 
cal constant than without (sufficiently large-scale modes 
will never enter the Hubble radius). The modes which 
enter near time '7//2, where r]f is the asymptotic final 
conformal time, will have a significantly delayed entry 
time, and therefore we expect the very largest scale ten- 
sor anisotropics to be somewhat reduced at the latest 
times. 



E. Optical depth 

One final line of sight effect on the CMB anisotropics 
that we will consider is the time dependent optical depth 
due to Thomson scattering. Looking back, it is the rapid 
increase in scattering near the time of recombination that 
makes it possible to speak of a "surface of last scatter- 
ing" , where the primary anisotropics are emitted. At 
later times, reionization results in a time dependent at- 
tenuation of the amplitude of the power spectrum, which 
it will be important to quantify. 

In our fiducial model the epoch of reionization occurs 
at redshift zji = 11.1, or an = 0.083, with the scale 
factor normalized to unity today. Here, we see a decrease 
in the amplitude of intermediate to small scale {£ ^ 30) 
anisotropics as photons are rescattered. On these scales, 
Ce [rj) is reduced by a factor e^"^" , where th is the optical 
depth to reionization. To compute the suppression in 
Ci{ri), recall the definition of the optical depth r(a. Gobs) 
between scale factors a and Oobs, given by 



(47) 



r(a,aobs) = ctt 



-; — nJa) da . 
da' ^ ^ 



(49) 
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where ctt is the Thomson scattering cross-section and 
ne(a) the electron number density. Assuming reioniza- 
tion is sharp and the energy density of radiation is sub- 
dominant, the optical depth to reionization, Tjf(aobs) = 
T{aR,aohs) is given by (see, e.g., [35]) 

rR(aobs) = 0.046(1 - Yp)-^ (JnJ^^a^Vw 

- ^fnJ^^a^^Vw^ , (50) 

with rfl(aobs) = for Oobs < Ofl- In this expression Yp is 
the primordial helium fraction (assumed to be 0.24). 

Eq. (50) gives Tfl;(l) = 0.088, so reionization reduces 
Ciiv) by ^ factor of approximately 0.84 by today [of 
course we must evaluate the spectrum at different times 
at the £ values related by the scaling relation Eq. (28)]. 
Much of the suppression of Ci{ri) occurs before the 
present time. Since the time the scale factor was half its 
present value, for example, the optical depth to reion- 
ization has only increased by 0.003. As Oobs oo, 
Eq. (50) says that the optical depth will only increase by 
Stoc = T]i{oo)—Tjf {l) = 5.7x 10^^, so that Ce{'i'i) will only 
be reduced by 0.1% relative to the value today. (Note 
that the peak positions should not be effected by reion- 
ization, since the peak structure is at scales smaller than 
those corresponding to the particle horizon at rescatter- 
ing.) Therefore the Universe is essentially transparent 
on cosmological scales today, and we have already seen 
nearly the maximum amount of reionization damping. 

F. Time evolution from CAMB 

In order to confirm and extend the analytic results of 
previous subsections we have modified the CAMB soft- 
ware [8] to compute the CMB power spectrum at differ- 
ent observational times. To do this, we simply modify all 
routines such that wc can evaluate the transfer functions 
T(fc,£, tls) at different 7'ls, using the set of best-fitting 
cosmological parameters as measured today. The changes 
required to CAMB are straightforward — for the most part 
all that is needed is to change the TAUO variable to trick 
the code into using a different observational time. The 
parameters we use are those of a standard six parame- 
ter spatially flat ACDM model, given by Qch'^ = 0.104, 
Qbh'^ = 0.0223, h = 0.734, ns = 0.951, As = 2.02 x 10-^ 
and zr ~ 11.1, where Vn{k) = As{k/ko)^^~^ with pivot 
scale ko ~ 0.05 Mpc""'^, and zr is the redshift of reioniza- 
tion. Where necessary, we set 715 = 1 to simplify compar- 
ison with analytic results. Wc do not expect significant 
differences in our results if these parameters are varied 
somewhat. Since the parameters are defined relative to 
their present day values, the spectrum measured by an 
observer with z > zr, for example, will not be affected 
by the epoch of reionization and a future observer will 
see a universe completely dominated by dark energy. 

In Fig. 2 we plot the power spectrum £{£+l)Cg{r]) / {2tt) 



calculated from our modified version of CAMB, parame- 
terizing the time dependence in terms of the observa- 
tional scale factor Oobs, where Oobs — 1 corresponds to 
today. It is clear that the angular scale of features re- 
sulting from projection of inhomogeneities near the LSS 
can be understood from the scaling relation (29). For 
example, in Fig. 2 we show the predicted scaling of the 
first acoustic peak position (located at £ = 221 at the 
present time) and find extremely good agreement with 
the predicted value. The existence of a future event hori- 
zon means that during A domination drLs / duohs tends to 
zero and so the acoustic peak positions become "frozen 
in". With our parameters we see that the first acous- 
tic peak becomes frozen in at the value £f — 290 as we 
predicted in Eq. (32). 




0.5 1 1.5 2 2.5 3 3.5 



log,,[l] 

FIG. 2: CMB power spectrum £{£ + l)Ce/{2TT) as a function of 
multipole £ and scale factor of observation flobs, as calculated 
from our modified version of CAMB for our fiducial ACDM 
model. We also plot the analytic scaling of the first acoustic 
peak (thick dashed line) predicted by Eq. (29). 

Recall that the scaling relation Eq. (29) predicts not 
only how the angular sizes of features in the spectrum 
scale with time, but also that the magnitude of the power, 
£{£ + l)Ci{ri), remains constant into the future at, e.g., 
any acoustic peak. For late times, Oobs ^ 0.3, Fig. 2 in- 
deed confirms this prediction. Our fiducial model under- 
went reionization at ur = 0.083, and we predicted in Sec- 
tion III E that as a result the power spectrum should be 
attenuated by approximately 16% on all but the largest 
scales. Again, this is visible in Fig. 2. Recall that we 
predicted a negligible 0.1% reduction in Ci{7]) between 
today and the distant future. 

Also visible in Fig. 2 is a substantial increase in power 
at the largest scales at late times due to the increas- 
ing ISW effect in our ACDM model, which we discussed 
in Section IIIC. Indeed, for Oobs 5.0 the quadrupole 
power actually exceeds the power at the first acoustic 
peak. The ISW contribution converges as the observation 
scale factor Oobs approaches infinity, since in the integral 
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in Eq. (42), we have g'{r]) — > and ij rjf as Oobs oo, 
where ??/ is finite. The asymptotic form of the power 
spectrum at late times is plotted in Fig. 3, together with 
the current spectrum. The dramatic increase in the ISW 
contribution, as weU as the shift in peak positions pre- 
dicted in Eq. (32), are clearly visible. Note that the ISW 
contribution to the dipole power, though not shown in 
Fig. 3, asymptotes to Cl^^ {r]f) ~ 5 x 10^^ for out model, 
according to CAMB. This corresponds to typical dipole 
ISW ampHtudes af^ ~ ^Cfw ^ 7 ^ IQ-^. Although 
such amplitudes represent a large increase over the cur- 
rent ISW dipole, they are still considerably below either 
the current total measured dipole, Oio = 2.5 x 10"'^ (with 
the polar axis aligned with the dipole) , or the amplitude 
of the galactic orbital dipole oscillation, aio = 1.5 x 10~^ 
according to Eq. (8). 




FIG. 3: Anisotropy power spectra for scalars (top pair of 
curves) and for tensors (bottom pair, arbitrary scale), calcu- 
lated using our modified version of CAMB for today (dashed 
curves) and for the asymptotic future (solid). 

Fig. 3 also presents the gravitational wave contribu- 
tion to the anisotropy spectra today and in the asymp- 
totic future. In Section HID we described the expected 
behaviour of tensor modes in the future, which entailed 
the same geometrical scaling of the small-scale cut-off in 
the spectrum, as well as a decrease in power at the very 
largest scales. Both of these features are visible in Fig. 3. 

As a check on our custom modifications to CAMB, 
we plot in Fig. 4 the power spectrum from CAMB for 
Cobs = 2.0, as well as the corresponding curve calculated 
from a power spectrum generated for today, Oobs = 1, 
and transformed to Oobs = 2.0 using the scaling relation 
Eq. (29). Additionally, the spectrum calculated from the 
scaling relation includes the increased ISW component 
calculated from the analytical approximation Eq. (42). 
To facilitate the use of this analytical expression, the 
spectral index was set to 715 = 1 for these calculations. 
Since the curves coincide at all but the largest scales, it 
is clear that the scaling relation has accurately captured 
the evolution of Ci{r]). However, the ISW contribution is 
substantially overestimated, indicating the limitations of 



the approximation Eq. (43) involved in deriving Eq. (42). 
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FIG. 4: Anisotropy power spectrum for a^hs ~ 2.0 calculated 
with our modified version of CAMB (solid line) and using the 
scaling relation Eq. (29) and analytical ISW approximation 
Eq. (42) (dashed). We set ns ~ 1 for these calculations. 

To make further contact with the analytic results in 
previous subsections, in Fig. 5 we plot the difference 
SCi = Ci{a'^^^) — Ci{aohs) calculated using our modified 
version of CAMB between the power spectrum today, at 
o-ohs = I7 and at a future time, when a^j^^ = 1 + Sa, for the 
cases Sa = 10~^, 0.001, and 0.01. These curves exhibit 
very accurately the scaling with Sa predicted in Eq. (33), 
when we recall that 6a = HqSti for small 5t]. Slight depar- 
tures from this simple scaling arc evident at the largest 
scales, where the spectra are nearly flat and hence their 
precise shape sensitively influences the location of zeros 
in SCi. In Fig. 5 we also plot the analytical result cal- 
culated from the power spectrum today using Eq. (45). 
Again we find excellent agreement with CAMB at all but 
the largest scales. We also find reasonable agreement at 
low £, indicating that our ISW approximation is quite 
good for small time increments from today. 



IV. THE DIFFERENCE MAP POWER 
SPECTRUM 

A. Analytical time evolution 

In the previous section we found a very simple scal- 
ing relation to describe the time dependence of the CMB 
power spectrum, which, together with the ISW effect, 
thoroughly describes the evolution of the spectrum. How- 
ever, if we are interested in the best way to observe 
evolution in the CMB we might expect that observing 
changes in the actual sky map, or the a^mS, should be far 
more promising than looking for changes in the heavily 
compressed Ci power spectrum. Intuitively, as the shell 
f = ?'LS of the LSS grows in size, we expect the finest 
structures to change first, then the larger ones. As we 
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FIG. 5: Absolute value of the difference in the CMB power 
spectra^(£+l)|(5Cf|/(27r) between Oobs ~ land a^bs = l + 5a. 
The solid curves were calculated from our modified version of 
CAMB, and from top to bottom denote Sa — 0.01, 0.001, and 



10" 



The dashed curve was calculated using the analytical 



expression Eq. (45) for Sa • 
calculations. 
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We used Tig = 1 for these 



shall see, the difference between two sky maps measured 
at different times does indeed encode much more infor- 
mation than the Ci spectra, namely the correlations be- 
tween the two maps, although perhaps counterintuitively 
the magnitude of a change Ci will dominate over the dif- 
ference map power spectrum for small time intervals. 



Definitions 



Since the variances Ce{r]) and Ct{r}') will in general differ, 

the quantity C^^ is not the best measure of correlations, 

and the spectrum Dj^ measures not only the loss of cor- 
relations but also the change in variance Cg {ij) . Therefore 
we may consider instead the modified difference map 

Safni = \ ^TT^O.tmi'q') - aern{l]), (56) 

V c^n'y ) 

which normalizes the modes at rj' to have the same vari- 
ance as those at rj. Then we find 



(57) 



where we have defined the normalized correlation func- 
tion by 



(58) 



This normalized function is useful in that we have Cj^ = 
1, 0, and —1 for perfect correlations, no correlations, 
and perfect anticorrelations, respectively. Similarly, the 
quantity 1 — C^^ measures the loss of correlations alone. 
However, the spectrum D^^^ will still be useful, since it 
measures both the loss of correlations and the change in 
variance, so it might be expected to be more sensitive to 
changes in the CMB than the quantity 1 — C^''' . Also, 

through the definition (51), the quantity D^'' is more 
directly tied to observations. 



Consider two measurements of the a^mS at times i] and 
r]' and define the difference map by 

5alm=a|,rn{'t^)-atm{^l)■ (51) 

Using Eq. (21), we can readily calculate the statistical 
properties of the difference map. We find 

((5a<?™(Sa|,„j,) Df' ^ti'^mm', (52) 

where we define the power spectrum of the difference 
map, , by 

of ^Ci{v) + CeX7f)-2Cf, (53) 

and 

^rTzik)T{k,e,r^s)T{K£,r'^s)- (54) 

Note that the quantity {5aim5a*^,^,) is diagonal in I and 
m, and that C^'' = Ci{r]). 

The quantity C^'' is an unequal time correlator, i.e. a 
correlation function that relates the anisotropics at time 
r] with those at time r/, through 

Re{aim{r])ai'rn'iv')) = C'f'' ^«'<5mm'- (55) 



2. Time evolution — flat sky approximation 

In analogy with Eq. (33) for SCg, we can write the 
spectrum of the difference map for small increments in 
time Srj = rj' — rj as 

Df^{^^)iS,?. (59) 

Therefore the difference of the power spectrum, SCi, 
dominates over the power spectrum of the difference map, 
D^^ , for small enough Sr], since SCi is only proportional 
to the first power of Srj. Of course for a particular St] 
we must calculate the coefficients of Srj and (Sri)'^ before 
we decide which method is more efficient if we are inter- 
ested in a detection. The details of instrumental noise 
are important and this is discussed fully in [14]. 

Beyond the [Sr])^ scaling, it is much more difficult to 
obtain the detailed evolution of D^^ than it was for SCi. 
Even when we consider only the Sachs- Wolfe plateau con- 
tribution, for which Ai{k) = —1/5 and A2{k) = A3(fc) = 
0, the Bcssel integrals involved in Eq. (54) cannot be an- 
alytically solved. In fact, this problem is related to a 
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divergence that can be illuminated if we employ the flat 
sky approximation described in the Appendix. 

Under that approximation, which is valid over small 
patches of sky and replaces the discrete indices £ and m 
with the continuous two-dimensional vector £, and the 
polar coordinate r with a Cartesian coordinate x par- 
allel with the line of sight, we can readily calculate the 
quantity on the left-hand side of Eq. (52). Using 



Sa{£) = a(£, Xls) - a{£, xls) 



(60) 



to define the difference map, where xls and the 
comoving distances to the LSS at times 77 and 77', respec- 
tively, and using Eq. (A. 6) for the anisotropics, we find 
that the result is not diagonal in £. Rather, it contains 
terms proportional to the Dirac functions 6^{a£ — £') and 
S^{a~^£ — £'), where we have defined 



Note that the integrand in Eq. (65), which is exact 
apart from the flat sky approximation, is bounded in 
magnitude by the integrand in Eq. (A. 9) for the power 
spectrum C{£, rj), as we vary Sxj^s- In place of Eq. (58), 
the normalized correlation function becomes in the flat 
sky approximation 



C"'" {£) 



^C{l,iT)a-^C{a£,rj')' 



(66) 



which is bounded by \C'^'^ {£)\ < 1. In the limit of short 
time interval, (5xls 0, we have C^^ {£) 1, corre- 
sponding to perfect correlation. We also have C^^ {£) 
as Sxhs — > 00 (recall, however, that in a ACDM uni- 
verse, only finite conformal time is available into the fu- 
ture). 



XLS 



(61) 



Indeed this is not surprizing: in the flat sky approxima- 
tion, an anisotropy on angular scale £ at ccls corresponds 
to a physical mode with comoving wavevector component 
'^/xhS orthogonal to the line of sight. But Eq. (10) tells 
us that such a mode should share correlations with the 
same physical scale at x^g , which corresponds to the an- 
gular scale £x'^g/xLs- Such off-diagonal correlations are 
completely suppressed in the full spherical expansion, as 
we found. 

The relevant quantity to calculate in the flat sky ap- 
proximation is instead the power in the difference map 
deflned by 

S^{£) = a{a£, x'^s) - a{£, xls)- (62) 
Again applying Eq. (A. 6) for a{£,XLs), we find 



Sa{£)Sa* (£'))= D''i"'{£) 6^{£-£'), 



(63) 



which is diagonal in £, with the power spectrum of the 
difference map given by 

£»''"'(£) EE C{£, if) + a-'^C{a£, i) - 2C""''(£). (64) 

Here C{£,r]) is the flat sky approximation to the 
anisotropy power spectrum, given by Eq. (A. 9), and 
(jvn (^f^ jg correlation function given by 



C"" (£) = — / dk, 



'-LS 



fc3 



(65) 

where Tik, k^) is the flat sky transfer function, Sx-^a, = 
^LS ~ ^LSi and kx is the component of the comoving 
wavevector parallel to the line of sight. Eqs. (62) to (65) 
are the flat sky analogues of Eq. (51) to (54), respec- 
tively. (The continuous argument {£) will always distin- 
guish quantities in the flat sky approximation from the 
corresponding exact quantities, which are labelled with 
the discrete indices fm-) 



3. Special cases 

Armed with the above flat sky approximation, we can 
now calculate the difference map power spectrum and 
correlation function in some special cases. First, con- 
sider the short time interval case, Sx^^ — *■ 0. Expanding 
Eq. (64) in powers of k^Sx^^ ~ ^<5a;Ls /a;LS , where £= \£\, 
we find 



D"" (£) 



' ,,vnmnKkxrki ^^^^ 



for the power spectrum of the difference map at low- 
est order in ^febs/a^LS- This expression exhibits pre- 
cisely the time interval dependence that we predicted 
in Eq. (59). (Note that in defining the difference map 
through Eq. (62), we have fixed the observed transverse 
wavevectors at both observation times, so the integrand 
in Eq. (67) is independent of ^77.) 

The integrand in Eq. (67) resembles closely that for 
the anisotropy power spectrum in Eq. (A. 9), but with an 
extra factor of fc^ in the numerator. In fact, using the 
relation 



we can easily rewrite Eq. (67) as 



\ilC^^'+^\£,ri)~fC{£,ri) 
XLS J ^ 



(68) 



(69) 



Here 6*'-"^+^-' (€, 77) is the anisotropy power spectrum cal- 
culated using a modified primordial power spectrum de- 
fined by 



(70) 



where ko = Iq/xls is the "pivot scale" used to define 
the primordial spectrum. (The result for D"^^ {£) is, of 
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course, independent of the pivot scale chosen.) For the 
special case of a power law primordial spectrum V-]z{k), 
with scalar spectral index ns, the modified spectrum 



r. 



("s+2) 

n 



(fc) has spectral index ns + 2; hence our choice 



of notation. Eq. (69) says that, for small time incre- 
ments, the shape of the power spectrum of the differ- 
ence map is determined entirely by the actual anisotropy 
spectrum "blue tilted", i.e. €^C(£, 77), together with the 
spectrum C^"^+^-'(£, 77) calculated from a blue-tilted pri- 
mordial spectrum, both evaluated at the same time rj. 
Therefore we expect that generically the shape of the 
difference map spectrum D'^'^ {£) will be roughly that of 
a strongly blue-tilted version of the anisotropy spectrum 
C{£, 77). The height of the spectrum of the difference map 
is determined by the ratio 6xi^s/xhs- 

Next, we can specialize to the case of the pure scale- 
invariant {ns = 1) Sachs- Wolfe plateau, which is charac- 
terized by r(fc, kx) = Ai{k) = const and Vnik) ~ const. 
Eq. (A. 9) gives in this case 



(71) 



in agreement with the standard Sachs- Wolfe result, to 
order l/£. The normalized correlation function is then 



C"'" {t) 



COs(fc^feLs) 



(72) 



In the short time interval limit, ^(5xls/2;ls ^ 1, 
Eq. (67) becomes for the Sachs- Wolfe plateau 



Sx 



LS 



a;LS 



(73) 



Note that this last integral is logarithmically divergent, 
but this is just an artifact of our assumption of a scale 
invariant spectrum to arbitrarily small scales [39] . Equiv- 
alently, Eq. (69) cannot be applied in this case, because 
the Sachs-Wolfe integral diverges for ns > 3. In real- 
ity, damping within the LSS imposes an effective cut-off, 
with essentially no structure at wavenumbers above some 
value /cinax [40]. Replacing the infinite limits with zbfcmaxj 
we can evaluate the integral in Eq. (73) with the result 
(valid for I/xls < fcmax) 



X-LS 



In- 



2/cvi 



1 



(74) 



This means that the contribution to the difference map 
power from the Sachs- Wolfe plateau is independent of 
£, apart from a logarithmic correction. This is the £- 
dependence we expect for the Sachs- Wolfe plateau for the 
anisotropy power spectrum Ci from a strongly blue tilted 
primordial spectrum, with scalar index ns = 3, as wc pre- 
dicted above based on Eq. (69). Comparing Eqs. (A. 9) 
and (67) for the power spectra of the anisotropics and of 
the difference map, and recalling the expression Eq. (A. 7) 
for the transfer function, we see that the "monopole" con- 
tribution to the spectrum £)'''' {£) (the part proportional 



to Ai) is proportional to the dipole contribution to the 
spectrum C{£,ri) (the part proportional to A2). 

Finally, we note that we can evaluate Eq. (65) for the 
correlation function analytically for all Sxi^s for the case 
of a delta-source in /c-spacc, Vnik) = Vn^{k — k). Such a 
source will be very helpful in understanding the temporal 
behaviour of the normalized correlation function C"''' {I) 
at late times. The result for such a source is 



cos kx{i)8xi,s if^<fca;LS, 
if^>fca;Ls, 



(75) 



where 



£2 



(76) 



is the line-of-sight component of the source mode k cor- 
responding to the observed scale t. This result tells us 
that the normalized correlation function is initially (at 
Srj = 0) unity, as expected, and subsequently oscillatory 
in 5rj, with positive correlations alternating with anticor- 
relations, and each scale t oscillating at a different rate. 
The largest angular scales (smallest t) reach anticorre- 
lation first, followed by smaller scales. The peak scale, 
(. = fcxLS: never becomes anticorrelated. This behaviour 
can be easily understood with the assistance of Fig. 6, 
by noting that at the peak £ scale we have kx[t) = 0, so 
that the modes k which contribute to the peak £ scale are 
parallel to the LSS and hence cannot produce anticorre- 
lations. As i decreases, kx{f) increases, i.e. k contains an 
increasing component parallel to the line of sight, so the 
first anticorrelations occur earlier and earlier. If we con- 
sider sources at different scales fc, Eq. (75) tells us that 
the first anticorrelations occur earlier for smaller scales 
(larger A:), as expected. 



4.. Origin of the difference map power 

As we mentioned above, the power spectrum of the dif- 
ference map, Dj^ , contains two distinct contributions: 
the loss of correlations and the change in variance be- 
tween the two times of observation. To make this explicit, 
and to determine which contribution is more important, 
we can use Eqs. (53) and (58) to write 



Ciiv) 



o 



' + 2 fl - CT' 



(77) 



1 / sce y 

4 \Ci{v)J 
Sr/ ^ 



(78) 
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FIG. 6: Two source modes, fci and fc2, with the same 
wavenumber k. Diagonal lines indicate troughs (solid) and 
crests (dashed). The observer is towards the bottom. The 
solid horizontal line indicates the position of the LSS at ini- 
tial time 77. The dotted horizontal lines indicate the position 
of the LSS at the first later time that produces perfect an- 
ticorrelations with the initial time, so that hot spots (solid 
circles) line up with cold spots (open circles) and vice versa. 
Mode fei, which corresponds to a smaller observed i, reaches 
anticorrelation before mode fe2 (since Sxi < 5x2)- A mode k 
parallel to the LSS would never reach anticorrelation, while a 
mode parallel to the line of sight would correspond to ^ = 
(in the flat sky approximation). 



difference map defined in Eq. (62) is closely related to 
the "normalized" difference map defined in Eq. (56). 

One further contribution to the difference map arises 
if we consider the absolute temperature anisotropies ST 
rather than the relative quantity dT/T, where T is 
the mean temperature. Recall from Section TUB that, 
if we consider the absolute spectrum T^C^ instead of 
the relative quantity C^, then the difference S{T'^Ci) = 
T'^{f]')Ce{i]') ~T^{r])Ci{r]) receives an extra contribution 
due to the expansion redshift. In that case we showed 
that the extra contribution is of the same order as the 
geometrical scaling part [recall Eq. (35)]. 

We can now repeat this calculation for the difference 
map power spectrum. The difference map in absolute 
temperature units is 



S{Tai,n) = T 



5t] 



{any 



-aen 



(82) 



at lowest order in 5ri/(aH)^^, where we have used T = 
—HT. Therefore the corresponding power spectrum be- 
comes 



The first line above is exact, while in the second we have 
dropped higher order terms in 



5C, 



Si] 
r}LS 



(79) 



[recall Eq. (34)] . With a calculation similar to that lead- 
ing to Eq. (67), it is straightforward to show that, for 

short time intervals (drj/rji^s <C 1), we have 1 — C^^ oc 
{St] / rji^s)'^ , so that the two terms in square brackets in 
Eq. (78) are of the same order in Srj/rji^s- 

The first term in square brackets in the expression (78) 
is due entirely to the change in variance 5Ci, while the 
second term is due solely to the loss of correlations be- 
tween a^miv) a£m('7') [recall Eq. (57)]. But from 
Eq. (69) we have 



Si] 



(80) 



This expression dominates the change in variance con- 
tribution to Eq. (78) by a factor Therefore, for all 
but the very largest angular scales (smallest i), the sec- 
ond term in the brackets in (78) must dominate over the 

first, and so the power spectrum D^^ is dominated by 
the loss in correlations. This can be confirmed by a direct 
computation in the flat sky approximation, which gives 



2C(^, 77) ( 1 - C""''(^) (^), 



(81) 



at lowest order in (577/?;ls- This means that the flat sky 
approximation to D^'' captures only the (dominant) con- 
tribution due to loss of correlations. This is not surpriz- 
ing: because of the scaling relation (29), the flat sky 



{S{Taem)SiTae,nr) = 



SCi 



Ci 



, (83) 



where we have used the expressions (52), (53), and 
(55). Next, retaining only terms to lowest order in 
Si]/{aH)^^ ~ Sri/i]i^g, and using Eq. (79), we have 



{6iTai,n)SiTai,nr) = 



D 



vn 



O 



Cf, 



(84) 

But then Eq. (80) tells us that the first term on the 
right-hand side of Eq. (84) dominates for all but the 
very largest angular scales [just as we argued above for 
Eq. (78)], and so 



{S{Tagjn)S{Tai,n)* 



(85) 



In other words, the part of the power spectrum for the ab- 
solute difference map 5{Tairn) which is due to the expan- 
sion redshift is subdominant. Thus, contrary to the case 
with SCi, it is irrelevant for the difference map whether 
we consider absolute or relative temperature differences 
(apart from on the very largest scales). 

In hindsight this result could have been anticipated 
directly from Eq. (82), since we expect that the change 
Sa£m, corresponding to the time interval ^77, should be 



Sai 



Srj 



lat 



(86) 



so that the first term on the right-hand side of Eq. (82), 
which is due to the expansion redshift, is subdominant 
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on all but the largest scales. Intuitively, the change in 
agm due to a change in observation time Sr] grows as the 
wavelength of the source modes decreases (for constant 
mode amplitude), since the corresponding increase in ra- 
dius of the LSS is a larger fraction of a shorter wavelength 
mode. On the other hand, the change in Taem due to the 
expansion rcdshift is independent of scale £. 

Similarly, the contribution to the difference map due 
to loss of correlations, which is described crudely by 
Eq. (86); is expected to dominate over the contribution 
due to changing variance Ct, which is roughly indepen- 
dent of i, as we showed rigorously above. 



B. Time evolution from CAMB 

1. Power spectrum and correlation function 

We have computed the correlation function C^^ from 

Eq. (54) and the difference map power spectrum D^^ 
from Eq. (53) numerically using our modified version of 
CAMB to extract T{k,£,ri,s) at different tls (as outlined 
in Section IIIF). using the cosmological parameters of 

our fiducial ACDM model. In Fig. 7 we display D^''' for 
the times rj and rj' corresponding to today, Oobs = 1, and 
future times when a'^y^^ = 1 + da, for the cases Sa — 10""*, 
0.001, 0.01, and 0.1. For small increments Sa these curves 
exhibit precisely the quadratic scaling D^^ cx {Sr])'^ that 

we predicted in Eq. (59), and the slope of D]^ for small 
£ matches our analytical prediction for the Sachs- Wolfe 
plateau, Eq. (74). For large increments Sa the difference 
map power spectrum approaches the sum of the individ- 
ual power spectra as the correlation function decays to 
zero, as we expect according to Eq. (53). Generally, these 
curves exhibit the heavily blue-tilted form we predicted 
in the previous subsection, due to the more rapid loss of 
correlations on smaller angular scales. 

Also shown in Fig. 7 is the curve calculated from 
the flat-sky analytical expression, Eq. (69), for the case 
Sa = 10~^. This curve coincides extremely well with 
the numerical result for £ > 20. The departures at large 
scales are due to two factors. First, the flat sky approx- 
imation is poor at those scales. Second, Eq. (69) was 
derived under the assumption that all anisotropies were 
primary, which is not the case for the ISW contribution. 

In Fig. 8 we plot the normalized correlation function 
Cg^ , calculated using our modified version of CAMB for 
our fiducial ACDM model, between the set of a^mS ob- 
served at 77 and 77'. Here, 77 corresponds to an observation 
of the CMB sky today at flobs = 1 and rj' to an obser- 
vation at a^i^g = 1 + (5a, where we illustrate the cases 
Sa = 0.001, 0.01, 0.03 and 0.1. For the smallest inter- 
val Sa, we find very strong correlation between the two 
sky maps, as expected. The correlations fall off as Sa 
increases, with the sky maps becoming somewhat anti- 
correlated for intermediate intervals before C^^ decays 



10" 
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FIG. 7: Difference map angular power spectrum D^'^ calcu- 
lated from our modified version of CAMB (solid lines) for the 
two times corresponding to aobs ~ 1 and a^bs = 1 + Sa, for 
the cases (top to bottom) Sa = 0.1, 0.01, 0.001, and 10"*. 
Also shown is the analytical result (dashed curve) calculated 
using Eq. (69) for Sa = 10"". 



to zero at the largest intervals. 

The general features of the correlation function can be 
understood by considering the detailed arguments pre- 
sented in the previous subsection. For Sa ^ 1 the in- 
crease in the LSS radius corresponding to the interval 
Sa is Sri^s = H^^^Sa. For the case of Sa — 0.01, using 
Ha = 73 kms"^Mpc"\ we find Sri^s = 40.9 Mpc, cor- 
responding to a comoving wavenumber k = 0.15 Mpc~^. 
This wavenumber is much larger than the wavenumber 
of the first acoustic peak, given by fc = tt/sls = 0.021 
Mpc~^, where sls — 150 Mpc is the sound horizon at 
last scattering. Hence, for this Sa, we are essentially 
sampling the same set of inhomogeneitics which give rise 
to the first acoustic peak at both times, and so we expect 
fluctuations to be correlated on these scales. Indeed, we 
see from Fig. 8 for Sa = 0.01 that C^''' ~ 0.9 for the 
first acoustic peak scale, £ ~ 220. Extending this argu- 
ment, we expect that C^'' 1 as £ for fixed Sa, 
as the largest scale (smallest k) features should be most 
correlated, and of course we similarly expect C^'' 1 
as (5a ^ for fixed £. 

The presence of anticorrelations was discussed in Sec- 
tion IV A, where we derived the behaviour of the cor- 
relation function in the flat sky approximation for the 
case of a delta-source at wavenumber k = k. The re- 
sult. Eq. (75), exhibited oscillating positive and negative 
correlations, with the first anticorrelations occuring ear- 
lier for smaller scales (larger fc), as we have confirmed 
here for a reahstic spectrum using CAMB. Eq. (75) also 
described anticorrelations occuring earlier for smaller £, 
with k fixed. This behaviour is not visible in the actual 
correlation function plotted in Fig. 8, since the real pri- 
mordial power spectrum is far from being a delta-source. 
If we consider a small subset of k modes, e.g. those corre- 
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FIG. 8: Normalized correlation function C^^ between the 
CMB sky observed today (a = 1) and at a = 1 + Sa, cal- 
culated with our modified version of CAMB, for 5a = 0.001 
(top panel, solid curve), 0.01 (top panel, dashed curve), 0.03 
(bottom panel, solid), and 0.1 (bottom panel, dashed). An- 
ticorrelations are seen to develop as 5a increases, before the 
correlation function decays to zero for large 5a. 



spending to the fourth acoustic peak scale, then some of 
those modes will be aligned nearly parallel to our line of 
sight and hence produce early anticorrelations at small £ 
for some 6a (recall Fig. 6). However, there are many more 
modes due to power at smaller k that are still tightly cor- 
related at the same Sa and hence result in C^^ ~ 1 for 
small i. 



2. Sky maps 

Assuming Gaussianity, generating a single realization 
of a set of aim^ usually involves drawing each ag^iv) ii^" 
dependently from a Gaussian distribution with variance 
Cg{ri). With the correlation function C^' , we have a 
measure of the degree of correlation between agmS at two 
different times. Hence, given a set of aimS at the first 
time, the variance of the distribution at each time, and 
the correlation between them, one can generate a real- 
ization of a second set of agmS at some later time. 

Formally, we draw the new set of agmS from the likcli- 



with 



27r|C|i/2 



C 



r<vv' r'^/v' 



(87) 



(88) 



where Xi is a random 2-vector containing each agm coef- 
ficient at 77 and ry', and we have relabeled the variance of 
the agm{i]) distribution at each time by C^^ and ^ . 

We illustrate the likelihood function in Fig. 9 for the 
a2m and a^m modes, where j] corresponds to an observa- 
tion at Oobs = 1 and rj' to Cobs = 2.0. The a2m coefficients 
are more tightly correlated than their a^m counterparts, 
since for such a large 6a the correlation rapidly falls off as 
£ increases. It is also noticeable that the contours of the 
likelihood are slightly elongated vertically, due to the in- 
creased variance of agm{rj') on large scales resulting from 
the increasing ISW effect. 
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FIG. 9: Distribution from which the a^mS are drawn. Here, 
77 corresponds to aobs = 1 and rj' to Oobs = 2.0, and contours 
show the 2(T error ellipse. The distribution of a2m is shown 
by the solid contour, which has a correlated set of variances 
(72(7?) and C2 (??'). The distribution of asm, shown by the 
dotted contour, has a smaller variance at both times and these 
are much less correlated. 

Therefore, our method of generating CMB sky maps 
involves firstly generating a random realization at some 
initial time, and then generating all subsequent realiza- 
tions by mapping the a^mS using the correlation function. 
For large 6a, where the a^mS are uncorrelatcd, we are es- 
sentially selecting a completely new set of coefficients. 
For small 6a, C^^ approaches unity and the atmS map 
trivially according to agmirj) a/rnij]')- some in- 
termediate intervals, anticorrelation favours a reversal of 
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sign of the a^mS, i.e. hot spots are mapped to cold spots 
and vice versa. 

We generate maps using the HEALPix code [41] with 
"side = 512, corresponding to a pixel resolution of 6.87 
arcmin. We present a series of these maps in Fig. 10, 
plotting the fractional temperature fluctuation 5Ti/T at 
each pixel i. For presentational clarity we show a patch of 
sky covering ^ 1000 square degrees, and use modes up to 
^max = 1000. We generate the first map at Oobs — 1, and 
show subsequent maps at Cobs = l+(5a, where 5a = 0.001, 
0.01, 0.1, and 1.0. We also show the difference map for 
each observation relative to Oobs = 1- We have checked 
that the power spectra reconstructed from our simulated 
sky maps agree with the intended spectra to within sam- 
ple variance. Note that for the sky maps in Fig. 10 we 
did not use the actual WMAP data for the present time; 
rather, we simply generated a random initial map accord- 
ing to the required Ci spectrum. 

Visually, the 5a = 0.001 map is extremely similar to 
the initial map. The variance of the map, given by 

\ / map ^ 

is over four orders of magnitude higher than the differ- 
ence map variance. For 5a = 0.01, the primary temper- 
ature fluctuations have a variance around two orders of 
magnitude more than the difference map, and changes in 
small scale structure (from the initial map) are clearly 
apparent. 

For 5a = 0.1 and 1.0, the variance of the difference 
is actually larger then the temperature fluctuations at 
that time, and acoustic scale structures are visible in the 
difference. This is understandable from our discussion of 
the correlation function — at these times the correlation 
on all but the very largest scales has dropped to zero, so 
that the variance of the difference approaches the sum of 
the initial and final map variances [recall Eq. (53)]. 

Finally, in Fig. 11 we present a simulated sky map 
for the asymptotic future. This map clearly differs from 
today's map, with the dramatic increase in large scale 
power due to the ISW effect readily apparent. For this 
map, we derived the initial a(„i coefficients up to £ = 20 
from the WMAP Internal Linear Combination map [17] 
[42] (for £ > 20 we generated random initial modes in- 
stead of using the real data, since the normalized cor- 
relation is negligible on those scales at these very late 
times). 

High resolution versions of these sky maps, 
together with animations illustrating the evo- 
lution of the CMB sky maps, Ci spectra, 
and correlation functions, are available at 
http://www.astro.ubc.ca/people/scott/future.html. 



V. DISCUSSION 

We have systematically described the temporal evolu- 
tion of the CMB, beginning with the mean temperature 
and dipole, and then moving to the anisotropy power 
spectrum. We found that the evolution of the spectrum 
is described at all but the largest angular scales by a 
simple scaling relation. At large scales the ISW contri- 
bution grows to dominate even the first acoustic peak at 
late times. The extra optical depth due to reionization 
is negligible into the future. 

We have introduced a correlation function between the 
CMB sky maps at different times which quantitatively 
encodes the intuitive notion that for small enough obser- 
vation time intervals 5r] and for source modes with small 
enough wavenumber fc, the anisotropies observed at the 
two times should be very similar. Closely related is the 
power spectrum of the difference map D^^ . We showed 
that the difference 5Ci scales like Srj for small intervals, 
while Z)^''' scales like {5'/])'^. The sensitivity of D]^ to 
changes in the sky maps is dominated by the loss of cor- 
relations at small angular scales, and the contributions 
from the change in variance Ce, as well as the change due 
to expansion redshift (if we consider absolute quantities) 
are subdominant. All of our numerical results were in- 
dependently confirmed analytically, and the validity of 
the necessary analytical approximations was elucidated 
by the numerics. 

The quantities we described in this work will be crucial 
to answering the question of the experimental detectabil- 
ity of a change in the CMB, or, more precisely, the ques- 
tion "how long must we wait to be able to confidently 
observe a change?" While the different time interval scal- 
ings we found for 5Ci and D'^'' might suggest that at- 
tempting to measure 5Cg would be much more favourable 
for small 5r], the situation is more subtle. In a separate 
paper [14] we quantify the detectability of changes in the 
CMB. 

On sufficiently long time scales the CMB evolution we 
have described will be obvious, and hence it is natural 
to ask what cosmological information a measurement of 
such changes might eventually provide to future cosmol- 
ogists. A measurement of the cooling rate of the mean 
temperature T would provide an independent and novel 
determination of the local Hubble rate Hq, through the 
relation T = —HT. As far as the primary tempera- 
ture anisotropies are concerned, according to Eq. (69) 

the shape of the difference map spectrum D^''' is deter- 
mined by the spectrum [and the primordial spectrum 
through the quantity C^"^^'^^f.,r])], so a measurement 

of the shape of D^'' should not provide much new infor- 
mation. With our assumption of a spatially flat geom- 
etry, the amplitude of D^^ is determined by the ratio 
5xi^s/xhs = <5t/(aoXLs). Thus a measurement of the am- 
plitude would directly fix the LSS radius xls and hence 
provide an independent constraint on the parameters ri„i, 
r^A, and Hq. (This determination of a;LS is replaced by a 
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FIG. 10: Simulated map realizations (top and left panels) and difference map (relative to Oobs ~ 1) (right panels) for 5a = 10~^ 
(middle panels, corresponding to a 13 Myr interval) and 10"'^ (bottom panels, 130 Myr). Note the vastly different power scales 
between the sky and difference maps. The maps presented here are for a patch of sky covering ~ 1000 square degrees. High 
resolution version available at http://www.astro.ubc.ca/people/scott/future.html. 
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FIG. 10: (Continued.) Simulated map realizations (left panels) and difference (right) for Sa = 0.1 (top) and 1.0 (bottom). High 
resolution version available at http://www.astro.ubc.ca/people/scott/future.html. 



determination of the angular diameter distance for spa- 
tially curved models.) The radius of the LSS is currently 
fixed, through observations of the accoustic angular scale, 
only up to the uncertainty in the matter content at last 
scattering. At very late times much additional informa- 
tion will of course become available as new modes become 
visible on the growing LSS. 

We have focusscd entirely on primordial anisotropics 
here. There are additional issues which arise when 
one considers secondary anisotropics, like gravitational 
lensing and Sunyaev-Zel'dovich effects, as well as time- 
dependent foregrounds of course. Such considerations 
depend much more heavily on the less well understood 
non-linear scales of structure, and so we leave this for 
others to pursue. We expect that there is plenty of time 
to pursue these ideas before any of these variations would 
be detectable. 

When this work was nearly complete, a related study 
appeared by Langc and Page [36]. Those authors calcu- 
lated future Ce spectra using CAMB. They also defined a 



correlation fuction equivalent to our Eq. (54), calculated 
it using CAMB, and generated simulated future sky maps. 
While they made no attempt at an analytical description 
of the evolution, their numerical results appear to agree 
with ours where they overlap. Furthermore, they made 
an estimate of the observability of the CMB evolution, 
which we have deferred to Ref. [14]. 
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APPENDIX: THE FLAT SKY APPROXIMATION 

The Bessel functions appearing in the various expres- 
sions relating primordial fluctuations to observed CMB 
anisotropics severely limit the extent to which analytical 
results can be obtained. However, a simple approxima- 
tion scheme, based on treating a small patch of the sky 
(and hence of the spherical LSS) as flat, allows us to use 
ordinary plane wave expansions and thereby to do "CMB 
without Bessel functions" . This small-angle approxima- 
tion is expected to be accurate up to terms of order 1/^, 
so that it is entirely appropriate for describing the acous- 
tic peak structure of the CMB. 

The flat sky approximation begins (see, e.g., [24]) 
by replacing Eq. (12) relating the observed temperature 
anisotropics with the perturbation functions on the LSS, 



in the strong coupling/free streaming approximation, 

by 



^ F{(l),{xLS,XLsd,VLs))- (A.l) 



Here is a 2-dimensional vector whose components rep- 
resent the angular displacement in two orthogonal direc- 
tions from the centre of the small patch of sky. In the 
Cartesian comoving coordinate vector (xLs,a:Ls^), the 
first component is parallel to, and the second two or- 
thogonal to, the line of sight to the centre of the patch. 
The coordinate value Xhs = V — '?ls refers to the comov- 
ing distance to the LSS from the point of observation. 
Analogously to Eq. (13) we can write 



d 

F{(l)i{xLs,XLsd, Vls)) = (^i(xLs, a^LS^, '7Ls) + (l)2{xLs,XLsd, Vls) (A. 2) 

for the monopole and dipole contributions. In place of the spherical harmonic expansion for the temperature fluctu- 
ation, Eq. (15), we here use a 2-dimensional Fourier expansion in terms of the continuous vector £ which replaces £ 
and m: 



T{r)) 2tt 



d'-laitie'^". (A.3) 



The statistical properties of the coefficients a{t.) can be determined in a manner completely analogous to that used 
for the spherical case in Section III A. Fourier expanding the perturbations (f)i according to 

</'.(xLS,a;LS0) = J d2fcxdfc,0,(fc)e*'=-■^-^^e''=^^^^ (A.4) 
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where and fc_L are Cartesian components of the wavevector k parallel and orthogonal to the line of sight, respectively, 
allows us to identify 

£ = XLsk±. (A.5) 

This tells us that £, the flat sky approximation to the spherical indices £ and m, is directly proportional to the 
component of the LSS fluctuation wavevector orthogonal to the line of sight, and that the relationship scales with 
the conformal time (or comoving distance) to the LSS, exactly as expected. Since fc_L is only a component of the 
wavevector k, Eq. (A.5) encodes the familiar fact that the mapping from fc to ^ is not one-to-one — rather, a range of 
fc's is mapped to a particular £ value. 
Using these expressions, we find 

1 f°° 

ai£,7^) ^ dfc,7^(fc,,.g/xLs)TFs(fc,A:.)e*'===^^^ (A.6) 

V27r.TLg J-oo 



where the fiat sky transfer function is 

Tps{k,k,) = Ai{k)+ikMk), (A.7) 

and the Ai are again defined by Eq. (11). Finally, using 
the statistical properties of TZ encoded in Eq. (10), the 
equal-time correlation function of a{£) becomes 



{a{e,v)a*{£',v))^C{i,ij)6\e~i'), 



(A.8) 



where 



7^.(fc)|T(M )l^ 



Again, we find that the coefficients a{£) for different £ 
modes arc uncorrelated. 

Notice that the time dependence of C{£, rf) is carried in 
the prefactor 1 /x^g as well as in the terms containing fc_|_ 



through Eq. (A.5) (if £ is held constant), whereas in the 
spherical case, Eq. (24), the Bessel functions carry the 
time dependence. Also, the complete absence of oscilla- 
tory functions in Eq. (A.9) means that it will be much 
easier to evaluate the CMB spectrum in the flat sky ap- 
proximation than in the spherical case, both analytically 
and numerically. 

In particular, we can easily apply Eq. (A.9) to rederive 
the scaling relation Eq. (29). Eq. (A.5) tells us that k± 
is invariant under the transformation a;LS 



'-LS 



and 



£ £' = feLs/a^LS- Therefore, Eq. (A.9) immediately 
implies that 



£'^C{£',ii)=£^C{£,r^), 



(A.IO) 



where £= \£\, regardless of the form of the transfer func- 
tions Ai(k) or of the primordial power spectrum V-jiik). 
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